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ABSTRACT 

We report the direct measurement of orbital period decay in the double neu- 
tron star pulsar system PSR B2127+11C in the globular cluster M15 at the rate 
of (—3.95 ± 0.13) x 10~ 12 , consistent with the prediction of general relativity at 
the ~ 3% level. We find the pulsar mass to be m p = (1.358 ± 0.010)M G and 
the companion mass m c = (1.354 ± O.O1O)M . We also report long-term pulse 
timing results for the pulsars PSR B2127+11A and PSR B2127+11B, including 
confirmation of the cluster proper motion. 

Subject headings: binaries: close — globular clusters: individual (M15) — gravita- 
tion—pulsars: individual(PSRB2127+HA, PSR B2127+11B, PSR B2127+11C) 



1. INTRODUCTION 

Pulsars in binary systems with neutron star companions provide the best available 
laboratories for testing theories of gravity. To date, two such systems have been used for 
such tests: PSR B1913+16 (Taylor & Weisberg 1982, 1989), and PSR B1534+12 (Stairs 
et al. 1998; Stairs et al. 2002). Both are consistent with Einstein's general relativity (GR). 

The globular cluster M15 (NGC 7078) contains 8 known radio pulsars, the brightest 
of which are PSR B2127+11A (hereafter M15A), a solitary pulsar with a 110.6ms spin 
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period; PSR B2127+11B (M15B), a solitary 56.1ms pulsar; and PSR B2127+11C (M15C), 
a 30.5 ms pulsar in a relativistic 8-hour orbit with another neutron star (Wolszczan et al. 
1989; Anderson 1992). The Keplerian orbital parameters of M15C are nearly identical to 
those of PSR B1913+16, though the former did not follow the standard high mass binary 
evolution (Anderson et al. 1990). With our data set spanning 12 years, M15C now provides 
a similar test of GR. 



2. OBSERVATIONS AND ANALYSIS 

We observed M15 with the 305 m Arecibo radio telescope from April 1989 to February 
2001, with a gap in observations between February 1994 and December 1998 roughly corre- 
sponding to a major upgrade of the telescope. All observations used the 430 MHz line feed, 
with 10 MHz of bandwidth centered on 430 MHz. 

Observations up to January 1994 were made with the Arecibo 3-level autocorrelation 
spectrometer (XCOR), which provided 128 lags in each of two circular polarizations and 
506.625 fis time resolution. The autocorrelation functions were transformed to provide 128 
frequency channels across the band, and these data were dedispersed at the appropriate 
dispersion measure (DM) (Anderson 1992) and folded synchronously with the pulse period 
for each pulsar. Observations were broken into sub-integrations of 10 minutes for M15A and 
M15C, and 20 minutes for the fainter M15B. 

Beginning in January 1999, we used the Caltech Baseband Recorder (CBR) for data 
acquisition. This backend sampled the telescope signal in quadrature with 2-bit resolution 
and wrote the raw voltage data to tape for off-line analysis with the Hewlett-Packard Ex- 
emplar machine at the Caltech Center for Advanced Computing Research (CACR). After 
unpacking the data and correcting for quantization effects (Jenet & Anderson 1998), we 
formed a virtual 32-channel filterbank in the topocentric reference frame with each channel 
coherently dedispersed at the dispersion measure (DM) of M15C (Hankins & Rickett 1975; 
Jenet et al. 1997). The coherent filterbank data were then dedispersed and folded for each 
pulsar as for the XCOR data. 

The folded pulse profiles were cross-correlated against a high signal-to-noise standard 
profile appropriate to the pulsar and backend (Fig. 1) to obtain an average pulse time 
of arrival (TOA) for each sub-integration, corrected to UTC(NIST). The standard pulsar 
timing package tempo 1 , along with the Jet Propulsion Laboratory's DE405 ephemeris, was 
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used for all timing analysis. TOA uncertainties estimated from the cross-correlation process 
were multiplied by a constant determined for each pulsar-instrument pair in order to obtain 
reduced x 2 — 1- An arbitrary offset was allowed between the XCOR and CBR data sets 
to account for differences in instrumental delays and standard profiles. The timing models 
resulting from our analysis are presented in Table 1, and post-fit TOA residuals relative to 
these models are shown in Figure 2. All stated and depicted uncertainties correspond to 
68% confidence. 



3. DISCUSSION 

3.1. Post-Keplerian Observables for M15C 

In addition to the five usual Keplerian orbital parameters, in the case of M15C we have 
measured three post-Keplerian (PK) parameters: advance of periastron (cu), time dilation 
and gravitational redshift (7), and orbital period derivative (P&). The dependence of these 
PK parameters on the Keplerian parameters and component masses depend on the theory of 
gravity; in GR, these relations are (see Taylor & Weisberg, 1982; Damour & Deruelle, 1986; 
and Damour & Taylor, 1992): 

/ p \ -5/3 

Cj = 3G 2 /VMgj {l-e 2 Y l M 2 /\ (1) 

7 = G 2 ^c- 2 e(^y 3 m c (m p + 2m c )M^, (2) 

P b = -±^G 5 /V 5 (g)~ 5/3 (l " e 2 )- 7/2 (l + ^e 2 + |e 4 ) m p m c M-V3 , (3) 

where G is the gravitational constant and c is the speed of light. The measurement of any 
two PK observables determines the component masses under the assumption that GR is the 
correct description of gravity; measuring the third parameter overdetermines the system and 
allows a consistency test of GR. 



3.2. Kinematic Effects on Pulse Timing 

The rate of change of orbital period that we observe in the M15C system, (P b ) ohs , is 
corrupted by kinematic effects that must be removed to determine the intrinsic rate, (P b ) mt . 
Following the discussion of Phinney (1992, 1993) regarding the parallel case of kinematic 
contributions to P, we have 
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(. \ km 
Ph\ i>n / , 8 — cos b cos / \ ad ai 
— = cos 6 cos + +- -, 4 

P b ) cR \ l + 5 2 - 25cosb cos I J c c' v; 

where v = 220 ± 20kms _1 is the Sun's galactic rotation velocity, R = 7.7 ± 0.7 kpc is 
the Sun's galactocentric distance, 5 = rf/i? 0) /•£ is the proper motion, d = 9.98 ± 0.47 kpc 
is the distance to the pulsar (McNamara et al. 2004), and a\ is the pulsar's line-of-sight 
acceleration within the cluster. The first term in equation (4) is due to the pulsar's galactic 
orbital motion, the second to the secular acceleration resulting from the pulsar's transverse 
velocity (Shklovskii 1970), and the third to the cluster's gravitational field. 

Acceleration within the cluster may well dominate the kinematic contribution to Pb, 
but ai is an odd function of the distance from the plane of the sky containing the cluster 
center to the pulsar, and since we do not know if M15C is in the nearer or further half of the 
cluster, we must use its expectation value, d\ = 0. Phinney (1993) calculates a maximum 
value of |a;| max /c = 6 x 10~ 18 s -1 , too small for the observed P to provide a useful constraint. 
However, the unknown a\ still dominates the uncertainty of (Pb) kin ; we take the median value 
of 0.71 |a«lmax as tne uncertainty in ai (Phinney 1992). Evaluating equation (4), the total 
kinematic contribution is 

( P b \ = (-0.0095 ± 0.12) x 10~ 12 , (5) 
and subtracting this contamination from (P?>) obs yields the intrinsic value 

lp b \ = (-3.95 ± 0.13) x 10~ 12 . (6) 

3.3. Component Masses of M15C and a Test of General Relativity 

Solving equations (1) and (2) given the measured values of cu, 7, Pb, and e gives m p = 
(1.358±O.O1O)M , m c = (1.354±O.O1O)M , and M = m p + m c = (2.71279 ±O.OOO13)M in 
the framework of GR (Fig. 3). This result is consistent with, and more precise than, previous 
mass measurements for the neutron stars in the M15C system (Prince et al. 1991; Anderson 
1992; Deich & Kulkarni 1996). We note that these masses are consistent with the masses 
of double neutron star binaries observed in the field (Thorsett & Chakrabarty 1999; Stairs 
2004). M15 is a metal-poor cluster with a mean metallicity [Fe/H] = -2.3 (Sneden et al. 
1991), suggesting that the mass of neutron stars is not a strong function of the metallicity 
of their progenitors. 
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Our determination of a third PK parameter gives a test of GR; (A) int is 1-003 ± 0.033 
times the predicted value. While M15C provides an impressive test of GR, it is less stringent 
than the 1% Co-^-P h test provided by PSR B1913+16 (Taylor & Weisberg 1989) and the 
0.5% u-'-f-s test provided by PSR B1534+12 (Stairs et al. 2002), where s = sini is the shape 
parameter determined through measurement of Shapiro delay. We note that the uncertainty 
in the intrinsic orbital period decay is due almost entirely to the kinematic contribution, so 
further observations will not significantly improve our determination of (Pb) mt or the quality 
of the test of GR it allows. 



3.4. Proper Motion of M15 

The proper motions resulting from our timing analysis give absolute transverse velocities 
for M15A and M15C several times greater than the cluster escape velocity. The measured 
proper motions for these two pulsars and M15B are shown in Figure 4, along with four 
published proper motion measurements for M15 based on optical astrometry. The pulsar 
proper motions are all consistent with each other; their average is \i a = (— 1.0±0.4) masyr -1 , 
fi$ = (—3.6 ± 0.8) mas yr^ 1 . This result is in good agreement with the cluster measurement 
of Cudworth & Hanson (1993). 



3.5. Intrinsic Spin Period Derivatives 

If we assume that GR provides the correct description of gravity, we can use (Pb) obs to 
determine the total kinematic correction to Pb and hence, to P for M15C. We find 

kin / . \ kin 

19 „-l 





= (-8± 17) x lO^s" 1 , (7) 

GR. \ / GR 

which corresponds to a\jc = (4 ± 17) x 10 _19 s _1 . We now apply this correction to the 
observed value of P and find the intrinsic value assuming GR, (-P)qr — (0.00501±0. 00005) x 
10~ 15 . This intrinsic spindown rate allows us to improve upon the previous estimate of 
the pulsar's characteristic age and magnetic field strength (Anderson 1992); we find r c = 
(0.097 ± 0.001) Gyr and B suii = (1.237 ± 0.006) x 10 10 G. 

Our timing models for M15A and M15C include P (Tab. 1) which is unlikely to be 
intrinsic to the pulsars. For M15A in the cluster core, Phinney (1993) estimates the kinematic 
contribution to be \di/c\ = P/P < 10 _26 s _1 (80% confidence). This is significantly larger 
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than the observed P/P 
jerk resulting from the c 



core, we measure 



3 x 10 8 s , so the observed P is consistent with the expected 
uster's mean field and nearby stars. For M15C, far from the cluster 
P/P ~ 10 _28 s _1 . We note that our measurement of P in M15C is not 
of high significance (~2<r), and may be an artifact of the systematic trends apparent in our 
timing data (Fig. 2). 
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Table 1. Pulsar Parameters for B2127+11A-C 



Paramctcr a 




Pulsar 






B2127+11A 


B2127+11B 


B2127+11C 


Right ascension, ctj2000 


21 h 29 m 58!2472(3) 


21 h 29 m 58 , !632(l) 


21 h 30 m 01?2042(l) 


Declination, 5j2ooo 


+12°10'0T.' 264(8) 


+12°10'00^'31(3) 


+12°10'38'.'209(4) 


Proper motion in o, /i a (mas yr _1 ) 


-0.26(76) 


1.7(33) 


-1.3(5) 


Proper motion in S, (is (mas yr _1 ) 


-4.4(15) 


-1.9(59) 


-3.3(10) 


Pulse period, P (ms) 


110.66470446904(5) 


0.05613303552473(9) 


30.52929614864(1) 


Reference epoch (MJD) 


50000.0 


50000.0 


50000.0 


Period derivative, P (10~ 15 ) 


-0.0210281(2) 


0.0095406(6) 


0.00498789(2) 


Period second derivative, P (10 -30 s _1 ). . . . 


32(5) 




-2.7(13) 


Dispersion measure, DM (pc cm -3 ) 


67.31 


67.69 


67.12 


Binary model 






DD 


Orbital period, P;, (d) 






0.33528204828(5) 


Projected scmimajor axis, asini (lt-s) 






2.51845(6) 


Orbital eccentricity, e 






0.681395(2) 


Longitude of periastron, ui (deg) 






345.3069(5) 


Time of periastron, To 






50000.0643452(3) 


Advance of periastron, ui (dcgyr~ 1 ) 






4.4644(1) 


Time dilation & gravitational rcdshift, 7 (s) 






0.00478(4) 


Orbital period derivative, (P 6 ) obs (10~ 12 ).. 






-3.96(5) 


Weighted RMS timing residual (/is) 


58.9 


103.5 


26.0 


Derived Parameters 


Orbital period derivative, (P b ) int (1(T 12 ) . . 






-3.95(13) b 


Pulsar mass, m p (Mq) 






1.358(10) 


Companion mass, m c (Mq) 






1.354(10) 


Total mass, M — m p + m c (Mq) 






2.71279(13) 


Galactic longitude, I (deg) 


65.01 


65.01 


65.03 


Galactic latitude, b (deg) 


-27.31 


-27.31 


-27.32 


Transverse velocity, v± (km s _1 ) c 


210(70) 


120(230) 


170(40) 


Intrinsic period derivative, P; nt (10 _15 ) b . . 






0.00501(5) 


Surface magnetic field, B sur f (xl0 1 0G) b .. 






1.237(6) 


Characteristic age, r c (Gyr) b 






0.097(1) 



a Figures in parenthesis are 1 a (68%) uncertainties in the last digit quoted. Uncertainties are calculated from twice 
the formal error produced by TEMPO. 

b Corrected for kinematic effects as described in §3.2 

c Bascd on the distance measurement of McNamara et al. (2004) 
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Fig. I. — Average pulse profiles for three pulsars in M15. Average pulse profiles of M15 
A (top row), M15B (middle row), and M15C (bottom row) as observed with XCOR (left 
column) and CBR (right column). The negative dips preceding and following the main pulse 
in the XCOR pulse profiles are an artifact of the correlator's coarse 3-level quantization; this 
distortion has been corrected in the CBR data (Jenet & Anderson 1998). 
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Fig. 2. — Timing residuals of M15A, M15B, and M15C. Residuals for each pulsar versus 
observation epoch are shown in the top three panels, with M15C residuals versus orbital 
phase in the bottom panel. Filled circles represent XCOR observations, with open squares 
indicating CBR data. 
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Fig. 3. — M15C mass-mass diagram. The constraints on the pulsar and companion masses 
in GR from the measured values of 7 (long dashed lines), cu (solid lines), and intrinsic P h 
(short dashed lines) are shown. The allowed region in mass-mass space at the intersection 
of these constraints is denoted by a heavy line segment near the center of the plot. The 
shaded region is excluded by Kepler's laws. The intersection of the constraints from the 
three post-Keplerian observables indicates that the behavior of this system is consistent 
with GR. 
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Fig. 4. — Proper motion of M15. Measured proper motions of M15A, M15B, and M15C 
in right ascension and declination, with rectangular regions indicating the published cluster 
proper motion measurements of Cudworth & Hanson (1993) (CH93), Geffert et al. (1993) 
(G+93), Scholz et al. (1996) (S+96), and Odenkirchen et al. (1997) (0+97). 



